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Abstract 

Chaos, oscillations, instabilities, intermittency represent only some nonlinear examples apparent 
in natural world. These phenomena appear in any field of study, and advances in complex and 
nonlinear dynamic techniques bring about opportunities to better understand animal signals. 

In this work we suggest an analysis method based on the characterization of the vocal fold 
dynamics by means of the nonlinear time series analysis, and by the computations of the parameters 
typical of chaotic oscillations: Attractor reconstruction. Spectrum of Lyapunov Exponents and 
Maximum Lyapunov Exponent was used to reconstruct the dynamic of the vocal folds. Identifying 
a sort of of vocal fingerprint can be useful in biodiversity monitoring and understanding the health 
status of a given animal. 

This method was applied to the vocalization of the Cervus Elaphus Corsicanus, the Sardinian 
Red Deer. 
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I. INTRODUCTION 



The physical and physiological mechanism of sound production are important to under- 
standing mammal vocalization which ranges from periodic vocal fold vibrations to com- 
pletely aperiodic vibration and atonal noise. Between these two extremes, a large amount of 
phenomena have been observed and reported : biphonation, cycles, subharmonic and 
chaotic behavior. These behaviors can be predicted by theoretical models. For example, 
the two mass model fthe most accepted for mammal apparatus of phonation) can exhibit 
irregular oscillations P, ^- The apparatus of phonation can be investigated through the char- 
acterization of the animal vocalization, where vocal nonlinearity can be used. According to 
Tokudaj^ the nonlinear analysis of human speech signal has been carried out extensively, 
while nonlinear characteristics for animal voice signals have not yet been investigated. 

Using the methods of nonlinear time series analysis we wished to understand the me- 
chanics of the vocal folds starting from the vocalization time series. The characterization of 
the vocal signal as a chaotic time series can give important information on the health status 
of the animal, since the oscillation modes are related to the status of the throat tissues and 
to the strength of the animal. Furthermore tissues shapes of vocal apparatus are different 
among the animals and the characterization of several chaotic signals can be used in the 
monitoring of biodiversity. 

The last remaining populations of a sub-specie of the red deer: the Sardinian deer ( Cervus 
Elaphus Corsicanus) are found in the well preserved evergreen forest of Monte Arcosu in 
Sardinia (a protected area owned by WWF Italy). The Cervus Elaphus is the largest and 
most phylogenetically advanced species of Cervus. Head and body length is 1.65-2. 65m, 
tail length is 0.11-0. 27m, height at the shoulder is 0.75-0. 15m, and weight is 75-340 kg. 
The largest and strongest male generally has the largest harem. In order to maintain this 
position of superiority he must constantly keep the distance with rival males by bellowing 
out, and chasing off potential rivals who come near his females. After vocalizing, the largest 
remaining males size each other up, and if antler and body size are comparable, they battle 
for the females. Their antlers lock and each male attempts to forcefully push the other 
away. The strongest and most powerful male wins and secures a harem (group) of females 
for mating. In this work an extensive characterization of the vocalization of Cervus Elaphus 
Corsicanus is presented by means of Lyapounov exponents of the chaotic oscillation evidence 
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of registered sounds. 
II. MATERIAL 

A number of different signals corresponding to different sound emissions were considered. 
Only clear and low noise sound emissions have been analyzed, in order to focus exclusively on 
meaningful vocalizations, and to avoid spurious effect. The vocalizations were recorded from 
adult males in their natural environment and digitized with a sampling frequency of 22050 
Hz. Fig^b) shows a small portion of the analyzed signal and in Fig^d), the spectrogram 
(512 points FFT) of the signal is shown. 

Discrete Fourier Transform (Fig^c)) was used to perform a preliminary spectral analysis 
on vocalization units. The presence of regions with high density of unresolved frequencies is 
a necessary, even if not sufficient, condition for the occurrence of chaotic dynamical regimes 
Non- linear dynamics analysis were, therefore, was limited to signal units characterized 
by broad-band features in the frequency domain. Results reported in the present work refer 
to a single signal 0.420s long. The time series examined consists of a 9455 points sampled 
at 22050Hz. 



III. COMPUTATIONAL METHODS 

The analysis of the time series was performed using the software package TISEAnQ] 
(Time SEries ANalysis) [3], valued as the most well known and robust algorithm set for 
nonlinear time series analysis. Typical steps are attractor reconstruction from time series and 
the characterization of the chaotic dynamic by means of Lyapunov exponents and maximum 
Lyapunov exponent (MLE). 



A. Attractor reconstruction 



The attractor of underlying dynamics has been reconstructed in phase space by applying 
the time delay vector method 

Starting from a time series s{t) = [si, . . . , s^] the system dynamic can be reconstructed 
using the delay theorem by Takens and Mane. The reconstructed trajectory X can be 
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expressed as a matrix where each row is a phase space vector: 



X — [Xi, X2, . . . , Xm]'^ 



(1) 



where Xi = [si, Si+T, • • • , Si-(DB-i)T] and M = N - (De - l)T. 

The matrix is characterized by two key parameters: The Embedding Dimension De and 
the Delay Time T. The embedding dimension is the minimum dimension at which the 
reconstructed attractor can be considered completely unfolded and there is no overlapping 
in the reconstructed trajectories. If the chosen dimension is lower than De the attractor 
is not completely unfolded and the underlying dynamics cannot be investigated. Higher 
dimension was not used due to the increase in computational effort. 

The algorithm used for the computation of De is the method of False Nearest 
Neighbors^. A false neighbor is a point of trajectory intersection in a poorly reconstructed 
attractor. As the dimension increases, the attractor is unfolded with greater fidelity, and 
the number of false neighbors decreases to zero. The first dimension with no overlapping 
points is De- 

The delay time T represents a measure of correlation existing between two consecutive 
components of D^-dimensional vectors used in the trajectory reconstruction. Following a 
commonly applied methodology, the time delay T is chosen in correspondence to the first 
minimum of the average mutual information function [jlfl] . 

B. Lyapunov exponents 

Chaotic systems display a sensitive dependence on initial conditions. Such a property 
deeply affects the time evolution of trajectories starting from infinitesimally close initial 
conditions, and Lyapunov exponents are a measure of this dependence. These characteris- 
tic exponents give a coordinate independent measure of the local stability properties of a 
trajectory. If the trajectory evolves in a A^-dimensional state space there are N exponents 
arranged in decreasing order, referred to as the Spectrum of Lyapunov Exponents (SLE): 



Conceptually these exponents are a generalizations of eigenvalues used to characterize dif- 
ferent types of equilibrium points. 



Ai > A2 > ■ ■ ■ > A, 



(2) 
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A trajectory is chaotic if there is at least one positive exponent, the value of this exponent, 
said the Maximum Lyapunov Exponent (MLE) gives a measure of the divergence rate of 
infinitesimally close trajectories and of the unpredictability of the system and gives a good 
characterization of the underlying dynamics. 

Starting from the reconstructed attractor X, it is possible to compute with the method 



of Sano and Sawada 
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the SLE consisting of exactly n = De exponents. This method 
is a qualitative one, and in presence of a positive exponents, Ai, a more accurate method is 
necessary for the computation. 

The method of Rosenstein-Kantz [l^l. Il^ is used to compute the MLE from the time 
series. This method measures in the reconstructed attractor the average divergence of two 
close trajectories in the time dj{i). This can be expressed as: 

d,{{) = C,e^^(^^*) (3) 

where Cj is the initial separation. By taking the logarithm of both sides we obtain: 

Indj(i) =lnCj + Ai(«At) (4) 

This is a set of approximately parallel lines (for j = 1,2, ... , M) each with a slope roughly 
proportional to Ai. The MLE is easily calculated using a least-squares fit to the average line 
defined by 

2/(z) = i^(lnd,(z)) (5) 

where (■) denotes the average over all values of j. Figure E^d) shows a typical plot of 
(lnc?j(i)): after a short transition there is a linear region that is used to extract the MLE. 



IV. RESULTS AND DISCUSSION 



The signal considered was characterized by highly complex patterns in which different 
transients with both periodic and apparently aperiodic features were identified. The ap- 
parently random behavior of the numerical series, easily detectable with a simple visual 
inspection of the sound pattern, was confirmed by the power spectrum and spectrogram. 
Three different regions were put into evidence: at low frequencies, between and 70 Hz, 
a first distribution of unresolved peaks is present, a sharp peak is also present at 450 Hz, 
while a broad band of frequencies, ranging between 850 and 1500 Hz, is easily detectable. 
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The chaotic characterization was performed calculating the embedding dimension D^; by 
the false nearest method and in FigI2fa) the result of the computation is shown. The figure 
reports the fraction of false neighbors with respect to the embedding dimension and a value 
of D^; = 4 was found. The delay Time was considered as the first minimum of the mutual 
information function, and the value T = 8 was found. 

Starting from the time series the attractor was reconstructed using the delay method, 
and in Figl^b) a three dimensional projection of the attractor is shown. The structure of 
the attractors, related to the chaotic oscillation of the vocal folds, demonstrated that the 
irregular behavior observed in the time series was not due to noise. 

In order to completely characterize the chaotic nature of the vocalization, the Spectrum of 
Lyapunov Exponents and the Maximum Lyapunov Exponent Ai were evaluated. In Fig|21^c) 
values of the four exponents are reported and the presence of a positive exponent was 
detected. The accurate value of the MLE was computed by the Rosenstein-Kantz method 
and a value of Ai = 0.48 was found by a linear regression of the curves in the region between 
and 20 iterations. 

The Kaplan- Yorke fractal dimension Dl of the attractor equal to Dl = 2.58, confirms 
the high dimensional fractal qualities of the strange attractor. 

V. CONCLUDING REMARKS 

The analysis method proposed in this letter was applied to the vocalization of an adult 
male of Cervus elaphus corsicanus and put in evidence the chaotic behavior of the irregolar 
oscillations in the signal considered. A full characterization by means of attractor reconstruc- 
tion, Spectrum of Lyapunov Exponents, and Maximum Lyapunov Exponent was performed. 
A positive value of MLE was found. Future work aimed at identifying different individuals 
through the discussed parameters, will consist in the analysis of other vocalizations looking 
for a vocal fingerprint that may be useful in biodiversity monitoring. 
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Figures 
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FIG. 1: (a) The analyzed signal, (b) A portion of the signal showing the irregular nature of the 
vocalization, (c) The power spectrum shows the typical Spectral contents of an irregular signal: 
a broadband and a continuous one. (d) The spectrogram of the signal showing a fundamental 
frequency of 70Hz and other frequencies until 2000Hz 
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FIG. 2: (a) Computation of the embedding dimension by the False Nearest Method apphed to 
the time series. The fraction of false neighbors decrease to zero at a reconstruction dimension 
De = 4. (b) The attractor reconstructed by the method of delays: This highly structured trajectory 
indicates that the signal is chaotic and that the irregular motion is not a noisy one. (c) The 
Spectrum of Lyapunov Exponent showing the presence of a positive Lyapunov exponent and three 
negative exponents, (d) Computation of the Maximum Lyapunov Exponent by the Rosenstein- 
Kantz algorithm. The value of Ai is obtained by a linear regression of the curves in the zone 
between and 20 iterations. The value Ai = 0.48 was found. 
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TABLE I: Results of the analysis performed on the vocalization signal. A positive Lyapunov 
exponent and the value of Kaplan- Yorke dimension indicates the chaotic nature of the signal. 



Parameter Value 

Delay Time T 8 

Embedding Dimension De 4 

Maximum Lyapunov Exponent Ai 0.48 

Kaplan- Yorke Dimension Dl 2.58 
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